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I. INTRODUCTION 

There has been much discussion [IJ [2] concerning the dif- 
ference in results between the KARMEN [3] and LSND [3] 
experiments regarding appearance of electron antineu- 
trinos, v e , from muon antineutrino, P M , sources. Ini- 
tially, the concern was that the mass-squared differ- 
ence, Am 2 , characterizing the oscillation scale does not 
match up with the differences observed in atmospheric [5 
and solar [6 neutrino oscillations (with the latter now 
both confirmed and superseded by the results of the 
KamLAND[7] experiment). This concern was based on 
the assumption that only the three known active neu- 
trino flavors participate in the oscillations, so that the 
third value of Am 2 is determined by the other two val- 
ues. 

More recently, it has been recognized that light ster- 
ile neutrinos may exist and participate in oscillation 
phenomena [5J O EHj- Because multiple cycles have not 
been explicitly observed, this raises a serious question 
regarding an assumption in the analyses of all oscillation 
experiments to date, namely that the oscillation scales 
are sufficiently separated so as not to influence the the 
values extracted using the functional relations in simple, 
two-channel mixing. We have previously shown [8] how 
a reduced rank see-saw pdj [12] which couples several dif- 
ferent oscillations, leads to more complex phenomena, as 
was long ago recognized by Fermi and Ulam[TB] whenever 
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more than two oscillators are coupled. 

If there are, in fact, multiple paths which contribute 
to electron neutrino appearance from a muon neutrino 
source, then the shortest baseline (largest Am 2 ) oscil- 
lation would appear as excursions from a rising longer 
baseline 8 . The other paths would have independent os- 
cillation parameters Amf i and Am^-, where, without loss 
of generality, we treat mass cigenstate " 1" is the domi- 
nant contributor to the initial flavor and " 2" as the dom- 
inant contributor to the final flavor. Then the oscillation 
from initial to final flavor can be represented as 

Pivp. — ► v e ) = A 2 sin 2 (Am 2 2 a;) 

+ B 2 sin 2 (Am 2 i a;) 

+ C 2 sin 2 (Am^a;) + • • • (1) 

where x — 1.27 L/E (m/MeV) and the dots indicate 
that, in principle, more than one intermediate channel 
i may contribute. The coefficients must all be posi- 
tive semi-definite to ensure positivity of the appearance 
probability. The additional terms produce the rising 
baseline. [5] so that the problem may be viewed as an 
oscillation of the usual two channel type, but occurring 
over a rising baseline. Of course, for x sufficiently small, 
all of the sin 2 arguments can be simultaneously expanded 
and the appearance probability develops purely quadrat- 
ically, viz. 

P{v^v e ) ~ [A 2 (Am 2 2 ) 2 + B 2 (Am 2 l ) 2 

+ a 2 (Am 2 i ) 2 + ...]x 2 (2) 

We provide some simple illustrations, which have the ad- 
vantage of being able to improve the compatibility be- 
tween the KARMEN and LSND results in a way that 
implicitly makes predictions for the eagerly awaited re- 
sults from the MiniBooNE[16 experiment. 
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II. EXAMPLE WITH CP CONSERVATION 



We suppose that an oscillation from to v e occurs with 
a value of Am 2 consistent with the LSND allowed range. 
However, as was shown possible in Ref. |5], we assume this 
occurs with coupling to other (here unspecified) channels 
that produces a rising baseline for the two channel oscil- 
lation. Thus, the probability for detecting a i> e of en- 
ergy £ at a distance L from a source, P{v^ — > v e ) is 
given by the explicit part of Eq.Q above, and here we 
choose some particular examples for the parameter val- 
ues (where we absorb the factor of 1.27 into coefficients 
so that the values in the arguments of the sine functions 
are 1.27 x Am 2 in units of eV 2 .): 



Appearance Probability vs. L/E 

(scaled to nominal LSND signal) 




- ■ Very High Am" 

- - Two Channel (High Am 2 ) 
~ Two Channel (Low Am") 

- Three Channel (a) 
■ Three Channel (b) 



s = L/E (ni/MeV) 



^2ch;High - 0-0045 sin 2 (0.8 L/E) (3) 

= 0.0600 sin 2 (0.2 L/E) (4) 

^tTchS - 0.005 sin 2 (0.7L/i?) 

+ 0.001 sin 2 (0.3 L/E) 

+ 0.0025sin 2 (0.4 L/E) (5) 

P^md = 0.004 sin 2 (0.7L/i?) 

+ 0.005 sin 2 (0.2 L/E) 

+ 0.002sin 2 (0.5 L/E) (6) 

for the appearance rates in the two-channel and multi- 
channel cases respectively. Although the three-channel 
mass relation is satisfied in this example, we empha- 
size that the intermediate channel need not be an ac- 
tive neutrino (if light sterile neutrinos exist [3]) and fur- 
thermore, this relation need not have been satisfied if 
two different intermediate channels make the dominant 
contributions [T] . For very large Am 2 , rapid oscillations 
will average to a constant appearance rate independent of 
L/E, which we use to set a normalization of 0.0026 con- 
sistent with the scale for the signal reported by LSND[4J. 

For the multichannel cases, Fig. [T] reprises the charac- 
ter of the result in Fig. (2) of Rcf.|8j in the usual L/E 
terms. Note that the appearance probabilities are virtu- 
ally indistinguishable in the low L/E region covered by 
the KARMEN and LSND experiments, although wide de- 
viations occur in the larger L/E region that MiniBooNE 
can address. It is also interesting to consider an E/L 
plot of the same function as in Eq.([T|), along with the 
corresponding distributions for simple two-channel fits 
to the LSND experiment at high and low values of Am 2 
as shown in Fig. [2j Again, for a very high value of 
Am 2 , the limited resolution results in an averaged, flat 
distribution, experimentally indistinguishable from one 
independent of E/L. 

Both figures include an indication of the range of L/E 
(or E/L) over which each of the KARMEN, LSND and 
MiniBooNE experiments are sensitive. 



FIG. 1: Two channel mixing v e appearance probabilities from 
£/ M for three values of Am 2 compared with the two example 
three-channel rates discussed in the text vs. L/E . 

Appearance Probability vs. E/L 



(scaled to nominal LSND signal ) 
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x = E/L (MeV/m) 



FIG. 2: Two channel mixing u e appearance probabilities from 
£/ M for three values of Am 2 compared with the two example 
three-channel rates discussed in the text vs. E/L . 



III. EXAMPLE WITH CP VIOLATION 

In the discussion so far, we have not allowed for CP vio- 
lation. If that occurs, additional terms [2] arise from the 
imaginary part of the same product of four mixing ma- 
trices that produces the positive coefficients above from 
the real part. These terms are of the form 

AP(^ -> v e ) = D sin(2Am 2 2 a;) + E sin(2Am 2 ,a;) 

+ F sm(2Am 2 , i x) H (7) 

Note that the constraint of positivity of the coefficients 
does not apply to these terms. In fact, for "maximal" 
CP violation, in the sense of the conventional angle 
5 = 7r/2, it is straightforward to demonstrate that, de- 
pending on which way one represents the mass differ- 
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ences, some of the coefficients above must be negative 
semi-definite. Specifically, in the Particle Data Group 
(PDG) formulation[14 for exactly three channels 



D 

A 2 
D 2 
C 2 



= -E = -F = Si 2 S 2 3Sl3C23Cl2C?3 



L 12 J 13' : '23 L -13 

2 2 2 2 
S 12 S 13 S 23 C 13 



2 2/2 
, 12 c 12 c 13l C 23 



2 2 
S 23 S 13 



(8) 



where s\2 — sin(^i 2 ), etc. as usual. Clearly, if any one 
of the conventional CKM/MNS[T5] angles vanishes (or 
reaches ir/2) then the CP- violating parts vanish as they 
must. Note also that the positivity of the sin 2 terms is 
not affected. 

It is the positivity of the appearance probability that re- 
quires the relation above between D, E and F - these 
terms all change sign for the CP-conjugate channel. 
Therefore, the sum of the coefficients of L/E for small 
L/E must actually vanish: 



D x [Am 2 12 - Am 2 u + Am 2 2l ] = 



(9) 



which is guaranteed by the relations among the three 
mass differences. 

We provide one example with CP violation of an oscil- 
lation that agrees with both LSND and KARMEN and 
predicts that the signal in MiniBooNE may be smaller 
than the largest value expected from the LSND results: 

P^*" e) = 0.0025 (sin 2 (1.0 L/E)) - 0.001 sin(2.0 L/E) 

+ 0.0005 (sin 2 (3.0 L/E)) - 0.001 sin(6.0 L/E) 

+ 0.001 (sin 2 (4.0 L/E)) + 0.001 sin(8.0 L/E) 

+ 0.01 (sin 2 (0.5 L/E)) (10) 

where coefficients of all additional terms are assumed to 
be negligibly small (and we have again absorbed the fac- 
tor of 1.27 into the numerical parameters). Of course, 
none of the Am 2 values matches with those inferred from 
other experiments that do observe flavor oscillations, so 
the scenario here is viable only in the case that this set 
of oscillations is proceeding through neutrino mass eigen- 
states that are dominantly sterile neutrino states, with 
small flavor components. 

For the opposite CP process, applicable to KARMEN 
and LSND, the contributions from the imaginary part of 
the product of the U matrices (the sin(2Am 2 a;) terms) 
change sign, so we have 



°^" / " ) = 0.0025 (sin 2 (1.0 L/E)) 

+ 0.0005 (sin 2 (3.0 L/E)) 

+ 0.001 (sin 2 (4.0 L/E)) - 

+ 0.01 (sin 2 (0.5 L/E)) 



I- 0.001 sin(2.0 L/E) 
I- 0.001 sin(6.0 L/E) 
0.001 sin(8.0 L/E) 
(11) 



This formula improves the agreement between KARMEN 
and LSND, slightly, over the lowest Am 2 fit to both, and 



Appearance Probability vs. L/E 

(scaled roughly to nominal v — >V c LSND signal ) 




FIG. 3: Four channel mixing v a appearance probability from 
i/fi for Eq.( |10[ ) and the appearance probability for the CP 
conjugate channel for Eq.(ll| as given in the text vs. L/E 



compared with 2 channel descriptions described previously. 



Appearance Probability vs. E/L 

(scaled roughly to nominal v — >v LSND signal) 



— ■ Very High Am" 

■ ■ Two Channel (High Am") 

— " Two Channel (Low Am") 
4 Channels + CPV 

Opposite CP: v — > v 




x = E/L (MeV/m) 



FIG. 4: Four channel mixing v a appearance probability from 
9ft for Eq.(10l and the appearance probability for the CP 
conjugate channel for Eq.(ll| as given in the text vs. E/L 



compared with 2 channel descriptions described previously. 



does so without requiring large or v e disappearance 



at larger L/E 
than 1%. 



the maximum loss is only slightly more 



We plot the formulae in Eqs.(10l and (111 vs. L/E 
and E/L in Figs. [3] and [4] for comparison with the CP- 
conserving oscillations shown earlier. As is apparent both 
from these Figures and from the crude Am 2 values, these 
formulae do not represent a "fit" to the data, but are 
simply an indication of the possibilities available if one 
does not arbitrarily constrain the entire neutrino oscilla- 
tion picture to three active Majorana neutrino flavor and 
(light) mass eigenstates. 
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IV. DISCUSSION 

The form in Eq.Q and the additional terms in Eq.Q 
have far too many parameters to be tightly fit with data 
available from present neutrino experiments, hence the 
predeliction for fitting to two-channel mixing scenarios. 
(Even when more channels arc attempted, the dominance 
of one scale has been assumed [5 or the two-channel fit 
results are used[T|.) However, since the rising baseline we 
observed possible [5] is roughly quadratic, corresponding 
to the opening of contributions from a longer wavelength 
oscillation, it should be viable to include the possibility of 
a rising baseline with one additional parameter, T, viz., 

Pipirtt -» v final ) ~ A 2 sin 2 (1.27 Am 2 L/E) 

+T '(1.27 Am 2 L/E) 2 + ■■■ (12) 

without CP violation, where T accounts for the rising 
baseline/additional channels, or even more compactly, 

P{i>i n it^Vfin) ~ S (1.27 Am 2 L/E) 2 

+V (1. 27 Am 2 L/E) 3 + ■ ■ ■ (13) 

where S absorbs all of relative amplitudes and ratios 
of Am 2 values of longer wavelength channels into one 
parameter and V absorbs all of the additional parame- 
ters for any number of CP-violating effects into another. 
In fact, these functional forms do not even depend on 
our initial physical ansatz , Eq.Q, and have only the 
disadvantage of not being applicable for large values of 
L/E. We recommend that all oscillation experiments test 
such functional forms to determine whether or not the \- 
squared per degree of freedom of the fit to their data is, 
or is not, improved by such additions to the standard 
two-channel analysis. 

As noted in Ref.[2], with two-channel mixing, the LSND 
and KARMEN experiments are in best agreement for low 
values of Am 2 because in that case, the difference in dis- 
tances affects the results quadratically in favor of LSND, 
whereas for very high values of Am 2 they should have 
seen the same size signal and hence arc in disagreement. 
However, here we see that, while even for an intermediate 
value of Am 2 the agreement would be marginal in two- 
channel mixing, the problem is reduced once the effect 
of a third channel on the baseline for the two-channel 
oscillation is included. The improvement is even more 
striking when CP violation is allowed. In fact, the ratio 
between the signal expected in the two experiments can 
achieve essentially the same value (or an even better one) 
as that obtained with a small Am 2 fit. 



The concern over a smaller value of Am 2 for LSND is 
the effect of the larger intrinsic mixing amplitude re- 
quired to match the data obtained in the region of small 
L/E: It predicts large effects, particularly disappearance 
rates, at much larger values of L/E typical of reactor 
experiments |1 7j . for instance. However, as our examples 
demonstrate, the rising baseline breaks the relation, seen 
in two-channel mixing, between the rate that an appear- 
ance signal increases with increasing L/E and the size of 
the signal in the initial range of the effect. As shown in 
our examples, the total appearance rate remains near 1% 
at all values of L/E, completely consistent with the lim- 
its from short baseline disappearance experiments 1 1 71 ITS] 
for both i/fj, and v e . We emphasize that this is true even 
though a two-channel fit would require a much larger 
overall amplitude in order for the signal to have grown 
to the size reported by LSND yet have remained too 
small to be observed by KARMEN with its shorter base- 
line (which one would think not likely to be significantly 
shorter). 

Finally, we note that the inclusion of CP violation, which 
is to be expected, (but not CPT violation, which would 
be revolutionary) further improves agreement between 
KARMEN and LSND and makes explicit testable pre- 
dictions for the results of MiniBooNE, as long as the 
possibility of light sterile neutrinos is allowed. We reiter- 
ate that small amplitude mixing to light sterile neutrinos 
poses no conflict with any known laboratory experimen- 
tal data. 

We conclude that only when oscillation experiments can 
all provide unbiased L/E distributions, rather than re- 
porting parameters for two-channel fits to the oscillations 
observed, will definitive conclusions be possible, regard- 
ing neutrino mass and mixing parameters, that are inde- 
pendent of theoretical biases. 

The examples we have presented also suggest that, at low 
energy, the MiniBooNE experiment may observe a con- 
siderably larger or smaller signal for v e appearance than 
would be expected from the two-channel fits to KAR- 
MEN and LSND. However, the examples also show that 
this would not necessarily contradict the results of either 
of those two experiments, whether considered separately 
or jointly. 

We thank Bill Louis and Geoff Mills for several discus- 
sions on this subject. 
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